toonl&s odd 

NEUROLOGY 



Potential use and challenges of functional connectivity 
mapping in intractable epilepsy 

Robert Todd Constable 1 - 2 - 3 - 4 *, Dustin Scheinost 1 - 3 , Emily S. Finn 4 , Xilin Shen 1 , Michelle Hampson 1 , 
F. Scott Winstanley 2 , Dennis D. Spencer 2 and Xenophon Papademetris 1,3 

' Department of Diagnostic Radiology, Yale School of Medicine, New Haven, CT USA 

2 Department of Neurosurgery, Yale School of Medicine, New Haven, CT USA 

3 Department of Biomedical Engineering, Yale University, New Haven, CT, USA 

4 Interdepartmental Neuroscience Program, Yale University, New Haven, CT, USA 



Edited by: 

Mark Holmes, University of 
Washington, USA 

Reviewed by: 

Kaspar Schindler, University of Bern, 
Switzerland 

Andreas Schulze-Bonhage, 
Universitatsklinikum Freiburg, 
Deutschland 

'Correspondence: 

Robert Todd Constable, School of 
Medicine, Yale University, 300 Cedar 
Street, The Anlyan Center, N132, New 
Haven, CT 06520, USA. 
e-mail: todd. constable@yale. edu 



This review focuses on the use of resting-state functional magnetic resonance imaging 
data to assess functional connectivity in the human brain and its application in intractable 
epilepsy. This approach has the potential to predict outcomes for a given surgical proce- 
dure based on the pre-surgical functional organization of the brain. Functional connectivity 
can also identify cortical regions that are organized differently in epilepsy patients either 
as a direct function of the disease or through indirect compensatory responses. Func- 
tional connectivity mapping may help identify epileptogenic tissue, whether this is a single 
focal location or a network of seizure-generating tissues. This review covers the basics of 
connectivity analysis and discusses particular issues associated with analyzing such data. 
These issues include how to define nodes, as well as differences between connectivity 
analyses of individual nodes, groups of nodes, and whole-brain assessment at the voxel 
level. The need for arbitrary thresholds in some connectivity analyses is discussed and a 
solution to this problem is reviewed. Overall, functional connectivity analysis is becoming 
an important tool for assessing functional brain organization in epilepsy. 
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INTRODUCTION 

Functional connectivity in human neuroscience refers to the syn- 
chrony of activity in anatomically distinct regions of the brain: if 
two areas are highly correlated in their activity over time, they are 
considered functionally connected. As measured by fMRI, func- 
tional connectivity relies on the blood oxygenation level dependent 
(BOLD) contrast mechanism (Ogawa et al., 1990, 1992), the same 
as that used in traditional task-based functional MRI studies. But, 
rather than examining changes in response to specific input stimuli 
in a block design or event-related paradigm, connectivity mapping 
can extract information from correlations in the fMRI time-course 
data while the subject is at rest, in the absence of any externally 
imposed task. Acquiring connectivity data in the resting state 
makes connectivity analysis easily adaptable to clinical scanning 
as it requires no subject participation other than to remain still in 
the magnet and is thus just like any clinical imaging study. It can 
also be used intraoperatively in anesthetized patients. 

Task-based fMRI provides exquisite maps of functional regions 
differentially involved in the execution of a specific task and is the 
foundation of most functional brain imaging research. Task-based 
fMRI as it is used in surgical planning in epilepsy is chiefly focused 
on identifying eloquent cortex that must be spared in a surgical 
intervention, while other techniques are used to identify the epilep- 
togenic tissue to be resected. Task-based fMRI is not well suited 
for identifying abnormal cortical or subcortical tissue throughout 
the brain, because only a very small number of regions typically 
show differential activation to any given task and the number of 
tasks that can be run in a reasonable time with sufficient statistical 



power is low. Task-based fMRI has not seen widespread clinical use 
(outside of mapping for surgical planning) because of limitations 
on the number of tasks, difficulties associated with presenting tasks 
and/or training subjects on the task in a busy clinical MR center, 
and the lack of whole-brain assessment from such studies. 

Functional connectivity data is different from task-based fMRI 
data in that it does not provide information as to which areas 
of the brain are differentially involved in the execution of a 
specific task, but instead provides a more basic measure reflect- 
ing how different brain areas are functionally connected to one 
another. As such, functional connectivity studies are potentially 
more appropriate for whole-brain surveys of functional abnor- 
malities such as the clinical challenge in epilepsy of identifying 
seizure-generating tissue elements or networks. The hope is that 
examining the network properties of the seizure-generating tis- 
sue will lead to a better understanding of the impact of epilepsy 
on the functional organization of the brain, assist in understand- 
ing comorbidities, and facilitate surgical planning, which in turn 
could lead to improved surgical outcome through identification 
and resection of the critical node(s) in the seizure-generating net- 
works^). Connectivity data could, for example, be used to identify 
brain areas to target with invasive intracranial recording elec- 
trodes. Resting-state fMRI is not limited to identifying tissue to 
be resected but also has the potential to predict cognitive change 
following resection if the resection plan involves specific nodes 
in functional networks. Research in all of these areas is ongoing 
and these topics are discussed in further detail in the sections that 
follow. 
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Functional connectivity mapping was first described by Biswal 
et al. (1995) in a seminal paper, revealing correlations between 
brain regions that are synchronized in time, through spontaneous 
fluctuations of activity, while the brain regions themselves may 
be quite distant. Brain regions with high resting-state temporal 
correlations are thought to be involved in the same network and 
the more well-formed the network, the stronger the connectiv- 
ity. Spontaneous fluctuations lead to changes in the BOLD signal 
and most of the work to date has focused on very low-frequency 
(<0.1 Hz) fluctuations. 

Spontaneous fluctuations have also been studied with both sur- 
face EEF and intracranial recordings (Towle et al, 1999; Andrzejak 
et al, 2001, 2006; Schevon et al, 2007; Ortega et al, 2008a,b; Lehn- 
ertz et al., 2011; Palmigiano et al., 2012) and local differences in 
the connectivity between nodes in the seizure focus area have been 
observed. For a review discussing both increases and decreases in 
synchronization and their role in epilepsy please see Jiruska et al. 
(2013). The EEC work, it should be noted, has primarily focused 
on higher frequency fluctuations in the 1^0 Hz range although 
several groups now are beginning to target much lower frequency 
oscillations. 

After the initial paper by Biswal et al. (1995), much effort was 
focused over the next decade on verifying that the correlations 
observed between brain regions indeed reflect true functional con- 
nections and not simply noise correlations due to physiological 
fluctuations. Without doubt, care must be taken to minimize the 
effects of physiologic noise as there are clear correlations asso- 
ciated simply with respiratory and cardiac signal fluctuations; 
however, the current consensus is that meaningful functional con- 
nections are found in continuously recorded BOLD fMRI data. 
Since about 2005, the field of functional connectivity mapping 
has been expanding rapidly through novel approaches to analysis, 
noise removal methodology, and applications to clinical popula- 
tions and basic neuroscience problems. In particular, the field is 
expanding to more in-depth analyses that move beyond simply 
examining correlations between two or a small handful of brain 
regions to capturing network information from a large array of 
nodes across the brain. 

For most functional connectivity studies, BOLD fMRI data is 
collected with the subject in the resting state with eyes either open 
or closed, and the patient is not required to perform a task. It is 
also possible to obtain connectivity data in the presence of a task 
and/or after brain-state manipulations, but it should be noted that 
the connectivity patterns can be slightly modified by task or brain 
state. Functional connectivity patterns have been shown to be sen- 
sitive to brain state as well as behavioral variables (Hampson et al., 
2006a,b; Johnson et al, 2006; Rogers et al, 2010; Bonelli et al, 
2012; Cole et al, 2012), and to vary with development (Fair et al., 
2008; Schafer et al., 2009; Myers et al, 2010) and age (Dosenbach 
et al, 2010; Hampson et al, 2012). 

Connectivity changes have also been reported in several clinical 
populations (Quigley et al., 2001; Lowe et al., 2002; Irwin et al., 
2004; Saini et al, 2004; Haas et al, 2006; Waites et al, 2006; Hoff- 
man et al, 2007; Schafer et al, 2009; Wang et al, 2009; Freilich and 
Gaillard, 2010; Myers et al, 2010; Bai et al, 201 1; Killory et al, 201 1; 
Zhang et al., 201 la,b; Bagshaw and Cavanna, 2012; de Groot et al., 
2012) and there are now more than 100 publications related to 



connectivity measures in epilepsy patients. There is evidence that 
correlations between time-varying BOLD signals reflect intrinsic 
functional connections in that they are present when subjects are 
both awake and under anesthesia (Vincent et al, 2007; Martuzzi 
et al., 2010, 201 1) and they are highly reproducible (Shehzad et al., 
2009). Overall, resting-state functional connectivity mapping has 
significant potential to reveal the functional organization of the 
brain and how it may be altered in different diseases or disorders. 

This review begins with an introduction to resting-state func- 
tional connectivity, including basic data collection and processing 
steps, the type of information that can be obtained, and various 
means of application. The particulars associated with resting- 
state methodology as applied in epilepsy are considered includ- 
ing strategies for voxel-based and region-of-interest (ROI) based 
analyses. Issues associated with selecting a connectivity thresh- 
old or avoiding such thresholds, ROIs, and the sensitivity of the 
method to choice of ROI are considered in addition to the emerg- 
ing field of connectivity-based parcelation for identifying minimal 
functional subunits for nodal analysis. Finally we provide a num- 
ber of early clinical results reflecting the potential of functional 
connectivity data to contribute to the clinical management of 
epilepsy. 

THE BASICS OF FUNCTIONAL CONNECTIVITY DATA 
COLLECTION 

Resting-state data is typically collected using a gradient-echo echo 
planar imaging (EPI) pulse sequence that provides whole-brain 
coverage with a temporal resolution of 3 s or less at a field strength 
of 3 T. Current state-of-the-art EPI involves the use of a multi- 
band/multi-plexed EPI sequences (Feinberg et al., 2010) that can 
provide whole-brain coverage with 2 mm 3 voxel dimensions and 
a repetition time (TR) of less than 1 s. Short TR in fMRI gen- 
erally provides better statistical power (Constable and Spencer, 
2001) and eases removal of physiological noise in connectivity 
data (Lowe et al, 1998). A minimum of approximately 5 min of 
such data is required to obtain connectivity maps with reason- 
able signal-to-noise ratios (SNRs). Such an acquisition easily fits 
within the constraints of a clinical MRI study, although in general 
acquiring more data is advantageous. 

Initial post-processing steps involve motion correcting the 
data (to ensure all volume acquisitions are aligned through 
time), and removing time-course signals that may be of no 
interest including the mean global signal through time, compo- 
nents associated with cerebrospinal fluid and white-matter signals, 
and the variables describing subject motion (the three trans- 
lation and three rotation directions). Temporal drift terms are 
often regressed from the data and a low pass filter (<0.1 Hz) is 
applied. 

The extent to which regressing the temporal fluctuations asso- 
ciated with motion is effective has been a topic of much discussion 
in the literature (Weissenbacher et al., 2009) and there is evidence 
for residual motion effects after correction (Power et al., 2012; 
Satterthwaite et al, 2012, 2013; Van Dijk et al, 2012). The issue of 
whether or not to remove the global signal mean also remains a 
hot topic as this can lead to the introduction of negative correla- 
tions in the data, making interpretation difficult (Fox et al., 2009; 
Murphy et al, 2009; Hampson et al, 2010). 
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Following the steps in the chosen preprocessing pipeline, 
functional connectivity is measured as the temporal correlation 
between the signal from any pair of voxels or ROIs. The following 
sections discuss the pros and cons of each approach in detail. 

CONNECTIVITY MEASURES 

Resting-state connectivity mapping has been used to examine 
functional connections between cortical regions since the first pre- 
sentation of the method by Biswal et al. (1995). In the sections 
that follow, we discuss how connectivity analyses can be applied 
to epilepsy for examining alterations in normal brain networks as 
a function of disease, or as a function of surgical intervention, as 
well as how connectivity mapping can be used to detect segments 
of tissue with abnormal functional connectivity with the goal of 
identifying the seizure-generating foci or network. 

ROI-TO-WHOLE-BRAIN 

Defining a ROI and performing ROI-to-whole-brain connectiv- 
ity analysis is probably the most common approach to exam- 
ining connectivity in the brain; this was the original method 
introduced by Biswal et al. (1995). Such an approach is moti- 
vated when an investigator is interested in a particular brain 
region (the seed ROI) and wishes to examine what other brain 
regions the seed is connected to, as well as how such connec- 
tions vary between healthy controls and a patient group, or pre- 
and post-surgical intervention (see Figure 1 for an example). 
In epilepsy this approach has been used to examine changes 
in language networks before and after anterior temporal lobe 
resection (Bonelli et al, 2012). In another study, Pereira et al. 
(2010) examined connectivity between the left and right hip- 
pocampi in mesial temporal sclerosis patients and found that 
relative to control subjects, patients with left hippocampal sclerosis 
showed a larger decrease in functional connectivity than patients 
with right hippocampal sclerosis in subjects with left-hemispheric 
language dominance. Morgan et al. (2012) performed a similar 
study using each of the left and right hippocampi as seeds for 
a seed-to-whole-brain analysis. This work revealed that the con- 
nectivity between the right hippocampus and the ventral lateral 
nucleus of the right thalamus could distinguish between seizure- 
free patients with left temporal lobe epilepsy (TLE) and right 
TLE patients and that in general, connectivity was greater in the 



seizure-free patient group with left TLE compared to the healthy 
controls. 

In another study of the medial temporal lobe, Pittau et al. 
(2012) selected four manually drawn ROIs and examined differ- 
ences in connectivity between 23 patients with mesial TLE and 
compared these to 23 age- and gender-matched controls in an 
ROI-to-whole-brain analysis. They found that patients with right 
MTLE had decreased connectivity between right amygdala and 
right hippocampus and the brain areas associated with the default 
mode network, some prefrontal regions, and contralateral mesial 
temporal structures. The left MTLE patients showed decreased 
connectivity between the left amygdala and left hippocampus 
to the default mode network, contralateral hippocampus, and 
bilateral limbic prefrontal regions. 

Moving outside the temporal lobe, a study by Killory et al. 
(2011) performed ROI-to-whole-brain analysis using eight ROIs: 
three defined used a functional task-based paradigm and another 
five spherical ROIs based on previous literature. They then com- 
pared the resulting maps from patients with childhood absence 
epilepsy to those of control subjects. An overall decrease in connec- 
tivity in the attention nodes was observed in the absence patients, 
who demonstrated impaired connectivity in the insula/frontal 
operculum and medial frontal nodes relative to healthy control 
subjects. 

In a study that used simultaneous EEG/fMRI to define the ini- 
tial seed regions, Negishi et al. (20 1 1 ) compared connectivity to the 
seed within the same hemisphere and the contralateral hemisphere 
to derive a connectivity laterality index. Results showed decreased 
laterality of functional connectivity in the patients that were 
seizure-free after surgery compared to those that had recurrent 
seizures after surgery. 

R0I-T0-R0I 

Another approach to connectivity analysis is to define a set of ROIs 
and to measure functional connectivity between all possible pairs 
of ROIs. This has been applied in numerous studies of epilepsy 
and some of these are summarized below. Once multiple ROIs are 
considered, network properties for these ROIs can be measured as 
discussed in the section on network theoretic measures below. 

In a study by Bettus et al. (2010), the investigators chose five 
ROIs in each hemisphere - primarily in medial temporal lobe 





FIGURE 1 | Seed-to-whole-brain connectivity mapping A reference BOLD 
signal time-course from a seed in Broca's area (top row) is correlated with the 
BOLD time-course for all gray-matter voxels in the image, revealing areas to 
which the seed is functionally connected (bottom row; hot colors = strong 
positive correlation, cool colors = weak or negative correlation). Such 



seed-to-whole-brain connectivity maps may be examined for a single epilepsy 
patient relative to a group of healthy control subjects, or for a group of 
patients with a similar pathology to a group of healthy control subjects, or 
correlated with behavioral measures. This approach and the results it 
produces are highly dependent upon how the initial seed region is chosen. 
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regions, defined using the Pick atlas tool in SPM (Maldjian et al., 
2003) - and examined connectivity between each pair of ROIs and 
compared this to data from 36 control subjects. They reported 
decreased connectivity in mesial TLE patients relative to con- 
trol subjects, with the largest decreases on the ipsilateral side and 
some increased connectivity on the contralateral side. In a similar 
study, Morgan et al. (201 1) defined hippocampal ROIs using both 
structural and functional information and revealed a relationship 
between functional connectivity and causal influence of the left 
and right hippocampi that varied with the duration of disease in 
a group of 19 mesial TLE patients. Looking at a broader network 
of nodes in a subsequent paper, Morgan et al. (2012) compared 
pre-operative fMRI connectivity data in two groups of patients 
who were postoperatively sorted into seizure-free and those with 
recurring seizures. They showed that the connectivity between the 
right hippocampus and the ventral lateral nucleus of the right 
thalamus could distinguish between seizure-free patients with left 
TLE from those with right TLE with high sensitivity and speci- 
ficity. The patients with recurring seizures after surgery generally 
exhibited different connectivity values in this network from those 
that were seizure-free. 

In a study of 11 children with intractable epilepsy, Widjaja 
et al. (2013) used independent components analysis to identify the 
default mode network and then calculated the functional connec- 
tivity between different pairs of ROIs in this network. They found 
reduced connectivity in the default mode network in children with 
medically refractory epilepsy. 

Using ROI-to-ROI connectivity analyses with a number of pre- 
defined ROIs, Bai et al. (201 1 ) investigated functional connectivity 
in a group of childhood absence epilepsy patients. In a unique 
twist on the ROI approach, they also examined ROI homologs in 
a cross-hemispheric study where the ROIs in this case were indi- 
vidual voxels. The primary finding was increased interhemispheric 
connectivity in the lateral orbito frontal cortex in the patient group 
relative to the healthy control subjects. 

A pattern classification approach was applied to another ROI- 
to-ROI study in epilepsy by Zhang et al. (201 la). This work exam- 
ined connectivity patterns between nodes in both hemispheres and 
found, much like the Negishi et al. (2011) study described in the 
previous section, that asymmetry in functional connectivity could 
correctly distinguish epilepsy patients from healthy controls with 
82.5% specificity and 85% sensitivity. 

NETWORK THEORETIC MEASURES 

While ROTto-whole-brain and ROI-to-ROI analyses can reveal 
much about functional connectivity, such approaches do not take 
the next step of considering nodes in the context of a functional 
network and the properties of that network as a whole. Recently 
there has been an explosion of interest in applying network theory 
(see Achard et al., 2006; Bullmore and Sporns, 2009; Bressler and 
Menon, 2010; Hagmann et al. , 20 1 0; Rubinov and Sporns, 20 1 0) to 
the analysis of functional connectivity data in order to characterize 
brain connections at both the nodal and network level. This allows 
for the observation of the networks associated with specific brain 
functions and generally moves fMRI from identification of indi- 
vidual nodes to systems involved in the execution of a task. Such 
theory, when applied to resting-state fMRI data, can characterize 



the topology of normal networks in the brain and, by extension, 
identify abnormal patterns of connectivity. 

Network theory measures have been applied in epilepsy to EEC 
or intracranial EEG data (Ponten et al, 2009, 2010; van Dellen 
et al., 2009; Varotto et al, 2012), MEG data (Chavez et al, 2010), 
and cortical thickness correlations (Bernhardt et al., 201 1), but to 
our knowledge only one published study to date has specifically 
examined network properties using fMRI data (Zhang et al., 201 1 a; 
although in the voxel-based approaches section below we include 
a small study by Stuffiebeam that uses the graph-theory mea- 
sure of degree to attempt to identify epileptogenic tissue). In this 
work, a total of 36 ROIs were defined using both functional (lan- 
guage and motor regions) and anatomic definitions [Brodmann's 
(1909) areas]. Five network measures - degree, strength, cluster- 
ing coefficient, closeness, and betweenness centrality (Wasserman 
and Faust, 1994) - were calculated for these nodes in individ- 
ual subjects as well as for input into a classification strategy to 
determine if such measures could distinguish medial TLE patients 
from healthy control subjects. The results confirmed that net- 
work measures such as these can indeed aid in classifying epilepsy 
patients and that the epilepsy process is associated with changes in 
network-level functional brain organization. 

Taken together, independent of whether the analysis is ROTto- 
whole-brain, ROTto-ROI, or at the level of network properties, 
these studies suggest that functional connectivity has potential 
for identifying disrupted circuitry as a function of disease and, 
perhaps more importantly, predicting outcomes from surgical 
intervention. In the sections that follow we outline some impor- 
tant issues encountered in these studies related to defining ROIs 
and determining thresholds for connectivity. 

R0I-BASED APPROACHES AND THE PROBLEM OF ROI 
DEFINITION 

As described above, many connectivity analyses require the a pri- 
ori definition of at least one ROI. While not often highlighted in 
previous publications, the choice of the seed ROI(s) and how the 
exact boundaries of that ROI are defined is critical. If an ROI con- 
tains multiple time-courses then the average time-course from the 
ROI may not properly represent any of the time-courses within an 
ROI and the results may be completely erroneous. Further, vary- 
ing the spatial definition of the seed can substantial changes in 
results. This is easily highlighted by observing that in a typical 
ROTto-whole-brain connectivity map (e.g., see Figure 1), there 
are often sharp transitions from positive to negative correlations; 
hence, moving the seed can result in a very different map. 

As is evident in the papers already discussed, numerous 
approaches to defining ROIs have been used to date. Task-based 
fMRI has been used to define specific functional circuits within 
which connectivity can be analyzed (Frings et al., 2009; Bonelli 
et al., 2012). This approach, however, suffers from the limita- 
tions of task-based fMRI studies in general in that only a very 
limited number of ROIs are activated by a task and thus whole- 
brain assessment of connectivity is not possible using such defini- 
tions. Another approach has been to use independent component 
analyses (ICA) (McKeown et al., 1998) to delineate brain regions 
(Luo et al, 2012; Mankinen et al., 2012) but these have typically 
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identified only a very limited number of networks - often, for 
example, fewer than 10. 

Anatomic ROI definitions have also been used extensively 
(Crespo-Facorro et al., 1999; Tzourio-Mazoyer et al., 2002; Makris 
et al, 2005; Shattuck et al, 2008; Zhang et al, 2011a). Such def- 
initions are ideal in structures that are well-defined anatomically 
(such as the hippocampus) but are difficult in areas such as the 
frontal and parietal cortices, and therefore the risk of mixing tem- 
poral signals into heterogeneous ROIs in these regions is high. 
Many investigators have used small spherical ROI placements 
(Shehzad et al, 2009; Bai et al, 2011; Bettus et al, 2011; Killory 
et al., 2011; Koyama et al., 2011); in this case the risk of mix- 
ing different functional time-courses decreases with the size of 
the defined sphere but is not eliminated. Many investigators have 
arbitrarily parcelated the cortex into anywhere from 100 to 1000 
nodes, but again, with such an approach, the node definitions may 
not necessarily reflect true functional boundaries. 

An emerging area of investigation involves performing whole- 
brain parcelation based on the time-courses themselves (van den 
Heuvel et al., 2008; Shen et al, 2010, 2013; Craddock et al, 2012). 
This approach appears very promising because it can provide min- 
imal functional subunits with uniform time-courses within each 
unit. An example obtained using the approach of Shen et al. (2013) 
is shown in Figure 2. Both Shen et al. (2013) and Craddock et al. 
(2012) have shown that ROIs extracted from these parcelations 
had higher functional homogeneity than anatomically defined 
ROIs and thus were more relevant for fMRI connectivity analyses. 
This parcelation approach using connectivity data itself appears 
to solve the problem of providing whole-brain ROI definitions 
for meaningful connectivity analysis. The next problem is how to 
apply such an approach to a patient population or to a group of 
patients. For example, if one generates a parcelation from healthy 
control subjects to investigate differences in network properties 
between control subjects and epilepsy patients, any results could 
be interpreted as due to actual differences in network properties, 



or a mismatch in the ideal functional boundaries for the parcela- 
tion nodes that come from imposing a control-derived parcelation 
on patients. This latter question can be addressed by directly com- 
paring a parcelation derived from the patient or patient group to 
the parcelation derived from the healthy control subjects. 

Thus, while it has been difficult to define ROIs for functional 
connectivity analysis, these new parcelation approaches appear to 
have solved this problem. Another approach, however, is to avoid 
the ROI problem completely by moving to voxel-level connectiv- 
ity analysis, in which each voxel in the gray matter is treated as an 
individual node and summary statistics on the network or con- 
nectivity properties of each voxel are obtained. This approach is 
described in the next section. 

VOXEL-BASED CONNECTIVITY ANALYSIS 

Voxel-level network analyses have been developed that can pro- 
vide insight into the functional connectivity of individual tissue 
elements, and a number of important studies have been published 
using such approaches (Buckner et al., 2009; Martuzzi et al., 201 1; 
Stufflebeam et al, 2011; Scheinost et al., 2012). Voxel-level analy- 
ses have the benefit of not requiring a priori definition of a ROI, 
instead treating each voxel as a node in a network analysis. In this 
approach a network measure such as degree can be calculated for 
each voxel. (Degree is the number of connections to a voxel above 
some arbitrary correlation threshold, i.e., r > t, where f = 0.25). 
Such a degree map, as shown in Figure 3 below, provides then 
for the first time a gray-scale contrast reflecting the functional 
connectivity of each tissue element and is a potentially power- 
ful approach to identifying regions that have abnormal functional 
connectivity on a whole-brain level. Such degree maps can also 
be used to define nodes for further ROI-to-ROI network analyses 
and/or to compare individual patients with control-group data. 

Such a voxel-based approach has been used to study brain 
changes in Alzheimer's disease (Buckner et al, 2009), to examine 
the effects of anesthetics on the human brain (Martuzzi et al, 





FIGURE 2 | A map of reproducible (across 79 subjects) functional 
subunits identified using resting-state connectivity data. These functional 
subunits are ideal nodes for connectivity analyses as they have highly uniform 



time-courses for each voxel within a given node by definition. This parcelation 
approach (Shen et al., 2010, 2013) provides a solution to the ROI definition 
problem in connectivity or network analysis of the brain. 
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FIGURE 3 | Gray-scale short axis MR images (from n = 42 healthy control 
subjects) with contrast reflecting the functional connectivity of each 
voxel as measured by the network measure of degree (brighter colors 



201 1) and to identify epileptogenic tissue in epilepsy (Stufflebeam 
etal.,2011). 

There are multiple approaches to calculating degree; three sim- 
ple schematics are shown in Figure 4 to illustrate the different 
approaches. In these examples the voxel-based degree measure 
can be calculated for connections encompassing the entire brain 
(whole-brain connectivity), within the same hemisphere as the 
voxel (ipsilateral connectivity) or spanning connections to the 
other hemisphere (contralateral connectivity). In each case the cal- 
culation is performed for each gray-matter voxel in the brain and 
the intensity of each voxel then reflects its number of connections 
above a predetermined correlation threshold. 

An example of such an approach applied in epilepsy is shown 
in Figure 5 below where the three degree measures are shown 
for a single patient relative to 20 healthy control subjects. The 
idea behind this approach is that seizure-generating tissue may 
have altered functional connectivity either because of the epilepsy 
processes themselves or because the tissue is not functioning nor- 
mally. Generally, in all three measures, we observe widespread 
decreases along with some increases in functional connectivity 
for the patient relative to the controls. Much needs to be learned 
about this approach, however, as there are subtle differences in 
the ipsilateral and contralateral measures that may provide rele- 
vant information for specific pathologies. In the case shown in 
Figure 5 the patient had suspected right medial TLE and indeed 



indicate higher degree). Such maps can be obtained for individual patients 
and compared to control-group data to isolate tissue elements with abnormal 
functional connectivity. 



there is a large region of decreased functional connectivity in 
the right temporal lobe relative to healthy control subjects across 
all three measures. However, there are many other regions that 
also show differences in functional connectivity relative to control 
subjects and it is currently an open question whether these differ- 
ences reflect part of the epilepsy network, brain reorganization of 
functional subunits as a result of having epilepsy, or some other 
neurophysiological mechanism. 

Given some clinical consensus via more conventional measures 
(surface and/or video EEG, other clinical data, and/or invasive 
recordings) there is often some indication that the patient has 
foci in a particular region or lobe. In the next example, shown 
in Figure 6 below, the right hippocampus was suspected to be 
part of the seizure-generating network and indeed the ipsilateral 
degree measure showed decreased functional connectivity in this 
region for the patient (see Figure 6A). (As an aside, we note that 
the ipsilateral degree measure is perhaps the most straightforward 
and easiest to interpret, whereas the whole-brain and contralateral 
measures are more difficult to understand in terms of the lateral- 
ization of the source of the problem.) As in the previous example, 
many other regions also show alterations in connectivity and some 
of these may be part of the seizure-generating network. Still, with 
a priori clinical information about the suspected site of seizure 
onset, it is reasonable to use the right hippocampus as a seed 
ROI for an ROI-to-whole-brain analysis. The nodes connected to 
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Whole Brain Ipsilateral Contralateral 

FIGURE 4 | Calculation of functional connectivity maps from of connections to that voxel above a correlation threshold across the 

resting-state fMRI data. The network measure of degree, reflecting whole-brain (bilateral; left panel), within the same hemisphere (ipsilateral; 

functional connectivity for a single voxel, is calculated by counting the number middle panel), or to the opposite hemisphere (contralateral; right panel). 





FIGURE 5 | Functional connectivity difference maps for a single epilepsy 
patient versus a group of healthy control subjects (red = connectivity 
increased in patient relative to controls; blue = connectivity decreased in 
patient relative to controls). Some of these regions are consistent with the 



seizure onset zones but clearly a large number of regions show altered 
connectivity. There are subtle differences between these three measures and 
more work is needed to determine the sources of these differences and the 
interactions with different pathological variants. 



this right hippocampal region are shown in Figure 6B with hot 
colors reflecting strong connectivity to the seed region. Combin- 
ing the ipsilateral connectivity map in Figure 6A with the map 
showing connections to the suspected focal region in Figure 6B 
using a logical AND operation, yields a much more circumscribed 
map Figure 6C, which highlights regions that have both altered 
functional connectivity and are part of the same network. (The 
most superior slices in the brain are not shown in the bottom 
two panels as there are no regions in these slices that satisfy both 
constraints.) 

Yet, despite their potential, voxel-based analyses of connectivity 
have not seen widespread application primarily for two reasons: 
(1) a lack of sensitivity, and (2) vulnerability to threshold effects. 

The first problem is partially solved by higher-field magnets 
of 3 T or more as well as the more recent move to ultrafast pulse 
sequences such as the multi-band/multi-plexed sequence devel- 
oped by Feinberg et al. (2010). In addition, sensitivity can be 



increased further by examining the distribution of degree values 
across the entire range of thresholds; this also addresses the second 
problem, as described in the next section. 

The second problem, the need for an arbitrary threshold to 
calculate the degree measure for each voxel, is unique to the appli- 
cation of network theory measures to fMRI data. In most network 
theory applications, such as tracking the transmission of a disease 
or analyzing friend links on Facebook, the decision tree is binary - 
either the disease is transferred or it is not; either two people are 
friends or they are not. However, in the case of fMRI connectivity, 
since connectivity is measured as a correlation, values may span 
a continuum from —1 to 1. Thus, an arbitrary threshold is typ- 
ically invoked to decide if two regions are connected or not. In 
general there is no principled way to choose this threshold and 
the results can change dramatically at different thresholds, which 
is a major issue for applying network or graph-theory measures 
(Rubinov and Sporns, 2010) to fMRI data. To get around this 
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FIGURE 7 | Regions of high intrinsic connectivity for left temporal lobe 
epilepsy (LTE) patients (n= 10) shown at different statistical levels At 

p < 0.05 corrected (top row) the ICD map (left) highlights significant patterns 
of increased connectivity in the left hippocampus whereas the degree map 
(right) does not show significantly different connectivity in this region. At a 
less stringent statistical threshold (bottom row), the degree map shows the 
same areas at p < 0.05 uncorrected (right) as the ICD maps (left) at either 
statistical threshold. This suggests that intrinsic connectivity is capable of 
detecting abnormal connectivity in regions of epileptogenic tissue and that 
the ICD approach yields higher sensitivity and reproducibility. 



FIGURE 6 | (A) Ipsilateral connectivity for a single patient versus a group of 
healthy control subjects showing multiple regions with altered connectivity 
(warm colors = higher connectivity for patient; cool colors = higher 
connectivity for controls). (B) In an individual analysis of the patient data, 
selecting the right hippocampus, which was the onset zone suspected 
from clinical data, as a seed region reveals a network of nodes with high vs. 
low connectivity to the right hippocampus. In (C) the intersection of the 
maps in (A) and (B) highlights regions that have abnormal connectivity 
compared to controls AND are part of the same network involving the right 
hippocampus. 



INTRINSIC CONNECTIVITY DISTRIBUTION ANALYSIS 

By characterizing the entire degree curve for any of the three voxel- 
based connectivity metrics described above (whole-brain, ipsilat- 
eral, and contralateral) , the need for a specific but arbitrary thresh- 
old is eliminated, providing both a more consistent measure of 
functional connectivity and a more sensitive measure of individual 
differences. Summarizing the entire degree distribution curve cap- 
tures information about all connections, weak to strong. This para- 
metric estimation approach therefore enables the interpretation of 
any differences (between an epilepsy patient and a control group, 
for example) to take the form of "more stronger connections," 
"more weaker connections," or "more even spread of connections." 
An example in a group of 10 TLE patients comparing degree at a 
single threshold (r > 0.25) and ICD is shown in Figure 7. 

The ICD approach improves on the earlier voxel-based degree 
measure through increased sensitivity as well as through more sta- 
ble interpretation of the results. To date no comprehensive study 
has been published using these emerging methods, but our lab as 
well as others are currently applying this to the problem of localiz- 
ing seizure foci or networks in epilepsy and no doubt more studies 
will appear shortly. 



arbitrary threshold choice, Scheinost et al. (2012) presented a new 
approach that examines the entire connectivity distribution curve 
across all correlation thresholds from 0 to 1. (The absolute value 
of the correlations can be taken, therefore taking negative corre- 
lations into account, or the positive and negative correlations can 
be treated in separate analyses). This new approach, referred to as 
the intrinsic connectivity distribution, is described in the section 
that follows. 



SUMMARY 

In summary, resting-state functional connectivity as measured by 
BOLD fMRI reflects intrinsic connections in the brain, provid- 
ing insight into how the brain is wired and how such wiring 
may be altered in disease or through surgical intervention. Like 
task-based fMRI, functional connectivity measures can be altered 
by task, drug or brain state, but unlike task-based fMRI, which 
reflects small changes in activity superimposed upon a very high 
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baseline activity level, functional connectivity measures intrinsic 
functional connections that exist even in the presence of deep 
anesthesia. 

The measures of intrinsic connectivity as described above, par- 
ticularly the voxel-based degree measure, provide for the first 
time a contrast mechanism in MRI based on function rather than 
anatomy. Since this approach requires no task, the data is acquired 
much in the same way that anatomic MRI data is acquired and thus 
these measures can be easily incorporated into clinical diagnostic 
radiology departments. 

Furthermore, while task-based fMRI provides functional infor- 
mation on only a very limited number of cortical regions 
(i.e., those few regions differentially activated by the task), the 



voxel-based connectivity measures described above provide func- 
tional information on the whole-brain without the need for a task. 
This makes it ideal for investigating the effects of pharmacologi- 
cal agents and the impact of diseases and specific pathologies on 
the brain without the need for a priori knowledge or selection 
of ROI. The use of functional connectivity methods in epilepsy 
will undoubtedly increase as we learn more about the functional 
changes that occur with epilepsy and/or as a function of surgical 
intervention. 
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